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Abstract 

We use spin correlations to constrain anomalous r-lepton couplings at the LHC including 
its anomalous magnetic moment, electric dipole moment and weak dipole moments. Single 
spin correlations are ideal to probe interference terms between the SM and new dipole- 
type couplings as they are not suppressed by the r-lepton mass. Double spin asymmetries 
give rise to T-odd correlations useful to probe CP violation purely within the new physics 
amplitudes, as their appearance from interference with the SM is suppressed by m,-. We 
compare our constraints to those obtained earlier on the basis of deviations from the 
Drell-Yan cross-section. 
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I. INTRODUCTION 


There exists a very active r-lepton program at the ATLAS and CMS collab¬ 
orations. TliCT have already studied the reconstruction of the Z-resonance in the 
di-tau mode [l|, l2[ and placed limits on physics beyond the standard model (BSM). 
Searches for heavy resonances such as Z' bosons have excluded such particles in the 
di-tau channel up to masses around 1 TeV There is also an active program 

to observe the Higgs boson in this channel [6|-l9(], as well as searches for additional 
neutral scalars that decay to In this paper we propose the use of spin 


correlations to constrain r anomalous couplings. 

An effective Lagrangian for BSM physics with a complete catalog of operators 

for the case in which the observed 
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up to dimension six exists in the literature 
126 GeV state is the Higgs boson in the SM. In a recent paper [l^ we studied the 
subset of dimension six operators that describe the r-lepton dipole-type couplings, 
as well as two dimension eight operators which couple tan leptons directly to gluons 
and are thus enhanced by the gluon luminosities [1^. In particular we discussed 


the constraints that can be imposed on these operators by studying deviations from 
the Drell-Yan cross section at the LHC as well as by bounding the cross section 
for production of r-leptons in association with a Higgs boson. We now extend that 
study to include spin correlations measured in the angular distributions of muons 
or electrons in leptonic decay modes. 

The couplings involved in the study are the r-lepton anomalous magnetic mo¬ 
ment and electric dipole moment given by al and d'l respectively, 


C = -T r 

and its corresponding weak dipole moments af and df, 


C = 


9 


2 cos 9w 


T (af iysdf) r 


( 1 ) 


( 2 ) 


The r-lepton dipole moments have been studied many times before in the literature 
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These anomalous couplings, Eq. [T] and Eq. |2l originate in the gauge invariant 


dimension six operators, in the notation of 11 


C 


■- 9 ^ <I>K. + 9 '^ 


'/iU 


h.c. 


( 3 ) 


The correspondence between these gauge invariant operators and the anomalous 
magnetic moment, electric dipole moment (EDM) and weak dipole moment (ZEDM) 
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of the leptons is given by^ 


x/2 


^ --Re {diB - d 


- a2 


IW) 


cb^ —-tt:—R e (dgyy + sin^ 9\y {d^B — diiY^'j 


j-y \/2 V , ■ y 

dl = ^ Im [diB — d^w) ) «£ = 


A2 

^/2 V 

A2 


Im {d(,w + sin^ 9w {deB — dew)) (4) 


where v ~ 246 GeV is the Higgs vacnnm expectation valne, 9w is the usnal weak 
mixing angle and A is the scale of new physics, which we take as 1 TeV for onr 
nnmerical stndy. 


As argned in [1^ the nsnal power connting for new physics operators is altered 
for dimension eight operators that conple a lepton pair directly to glnons dne to 
the larger parton Inminosities. This motivates onr inclusion of the “lepton-gluonic” 
couplings for the r in this study 


c = ^ [d^G + h. c. 


( 5 ) 


Here G^^ is the gluon held strength tensor and G^^'^ = its dual. 

If we allow for GP violating phases in the coefficients, drc and d^Q, the resulting 
gluon-lepton couplings take the form 


^ (R<i(<ixG) + Re«e) G-*'“'Gy 


\/2A4 


+ (lm(<ixG) G-^'-Gy + Im«G) G-^'-G^^) hd 


( 6 ) 


In Table H we summarize the Icr constraints that we obtained on the r-lepton 
anomalous magnetic moment, electric dipole moment and weak dipole moments 


assuming a 14% measurement of the Drell-Yan cross-section at LHC14 in Ref. 13 


We compare them to the best existing constraints from Delphi 27|, Belle 28| and 


Aleph 29|]. The results can be interpreted as a sensitivity to a NP scale A ~ 0.5 TeV. 


For comparison, the same measurement of the Drell-Yan cross-section constrains the 
NP scale of the dimension 8 gluonic couplings A ~ 1 TeV. The cross-sections we used 



mrO^ LHC-14 

niraX existing 

m-rdX LHC-14 

nirdY existing 

V = 7 

(-0.0054,0.0060) 

(-0.026,0.007) Delphi 

(-0.0057,0.0057) 

(-0.002,0.0041) Belle 

V = z 

(-0.0018,0.0020) 

(-0.0016,0.0016) Aleph 

(-0.0017,0.0017) 

(-0.00067,0.00067) Aleph 


TABLE 1: Summary of constraints for Icr bounds that can be placed on the r-lepton 
anomalous magnetic moment, electric dipole moment and weak dipole moments at LHC14 
from Ref. 


13l | compared to existing bounds. 


^ There is a typo in Eq. 4 of Ref. 131 that propagates to the conversion of our bounds from the 
gauge invariant basis to the tan anomalous couplings corrected in an errata. We also use here a 
different, more convenient, normalization for dew- 
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to obtain the constraints in Table [H are approximately qnadratic in the anomalous 
couplings indicating that the interference with the SM is very small. This is, of 
course, due to the fact that the interference between the SM and the dipole-type 
couplings is suppressed by the r-mass. 

In this paper we extend our previous study considering constraints that arise 
from spin correlations. These spin correlations evade the helicity suppression of 
the interference terms in the cross-section and produce observables linear in the new 
physics couplings. In this way it is possible to improve the constraints on the electric 
dipole moments and to study their CP violating nature through T-odd asymmetries. 


II. SPIN CORRELATIONS 


Spin correlations in r-pair production including anomalous dipole type couplings 
have been studied in Ref. |2lJ . In that paper, the spin density matrix for production 
of r-pairs in e’''e“ colliders in the CM frame was constructed and combined with 
the decay matrix for polarized r in its rest frame. That formalism exhibits the spin 
correlations explicitly but is not suited for our calculation. We want to construct 
(Lorentz scalar) correlations in terms of observable momenta at the LHC, namely, 
the muon (or electron) momenta and the beam momentum. Furthermore, we want 


to measure the correlations with event simulations using MadGraphS 3^. The main 


advantage of this approach is the ease in introducing different types of new physics 
with the aid of FeynRules j3l|. In this paper we limit ourselves to dilepton decays 
of the r pairs, but in a future publication we will address the hadronic decay modes. 


A. CP violating couplings 


The imaginary part of the effective couplings gives rise to electric and weak 
dipole moments of the r-lepton. These dipole moments are known to produce a 
double spin correlation linear in the anomalous coupling, of the form 


O 


2s 


m. 




U U Q. & 


( 7 ) 


This correlation originates in the interference between the CP violating edm ampli¬ 
tude and the CP conserving SM amplitude. In this case, however, the interference 
requires a fermion helicity flip and is therefore proportional to the r-lepton mass, 
resulting in a large suppression at the LHC. On the other hand, contributions that 
are quadratic in the new physics couplings do not suffer from this suppression, and 
Eq. [7] is useful to probe terms of the form ~ Re{drv) Im(d.rv), that is. 


O 


2s 




Z.'y u u ol & 


( 8 ) 


and similar terms proportional to both the real and imaginary parts of the couplings 
drG and d^Q. 
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With the muons (or electrons) in leptonic tan decay acting as spin analyzers 
this is measurable as 


O 


SS 


Pr+Pr-P^+pf- 


( 9 ) 


which requires at least partial reconstruction of one r-momentum direction and may 
be better suited for hadronic decay channels. 

To probe the anomalous couplings that violate CP with terms in the differential 
cross-section that are linear in lm{drv), but not proportional to the r mass, we resort 
to single spin correlations. For example, for the parton level process qq —)■ one 

hnds that the Z exchange diagram leads to the CP-odd correlation 


Ois ~ d^gA{i 


u) - P^TPt+Pt- A- 


— St+ 


( 10 ) 


where t, u and s are the parton level Mandelstam variables, qa is the axial vector 
coupling of the Z to the charged leptons and we have neglected the smaller vector 
coupling, Qv- Note that the parton momenta pi ^2 appear in a symmetric combination 
(from the two antisymmetric factors, {t — u), and the explicit {pi —P 2 )) and therefore 
this correlation does not vanish after the symmetrization of pi ^2 that follows from 
the convolution with the parton distribution functions for the LHC pp initial state. 

In order to write T-odd correlations that are sensitive to Eq. [inland are expressed 
only in terms of observable momenta we note that: 


• In leptonic r decay, the spin is analyzed by the muon (or the electron) mo¬ 
mentum. The simplest way to compute this is using the method of Ref. 
which shows that for leptonic r decay, Eq. [10] becomes 


O 


Is 


{i - u) 


Pt- ■ Pf,-e^^y^a,^iPlP2PT-Pl+ +Pr+ ■ P^,+ ^^^,u,a,pPlP2PAPT+ 


( 11 ) 


• In the lab frame at the LHC the r-leptons are highly boosted so their three- 
momenta are very close to that of the muons. Further, in leptonic r decay it is 
not possible to reconstruct the r momentum completely. We then replace the 
r momenta with the corresponding muon momenta in the lab frame obtaining 

0[i -^OC {i - u)^f,,u,aMP2PAP^+ (12) 


• In the lab frame, the sum and difference of the proton momenta are just the 
center of mass energy and the beam direction, and the two parton momenta 
appearing in Eq. [T2] have to be expressed in terms of these: 


Pi = 


P = Pl+P2 = 




( 1 , 0 , 0 , 1 ) 


\/5(l, 0,0,0) 


P2 = ^(1,0, 0,-1) 

qbe.m = Pi-P2 = \/5(0,0,0,1) (13) 


• This leaves us with two possibilities: 

[dbeam ' {Pfi+ Pti~^ Qbeam ' (P/x+ ^ 

^2 \Qbeam ' i.Pfj,+ T ) Qbeam ' ^ ^A‘~)]za6 (1^) 
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The CP properties of these two forms merit discussion. Since these correlations 
involve the beam momentum, a property of the initial state, they do not have definite 
transformation properties under CP. In fact, CP transforms LHC correlations 
into anti-LHC {pp collider) ones. If we consider the same correlations for a pp 
collider instead, we see that in this case the hrst one is CP-odd and can only be 
produced by the electric dipole moments. The second one, however, is CP-even 
and cannot be produced by electric dipole moments at a pp collider. This is a novel 
feature for colliders, that does not occur at the parton level where the pre-factor 
q ■ {Pt+ +Pt-) = (Pi — P 2 ) • (Pi +P2) vanishes. We illustrate this with an example in 
Table IVinl 

Noting again that the symmetry of the initial pp state at the LHC forbids terms 
linear in %eann we use the correlation 

Ctest — [qbeam ' {Pfi+ ^ PtJ‘~)]lab 

to gauge the statistical signihcance of our asymmetries. Interestingly, as seen in 
Table IVHIl this asymmetry does not vanish in pp colliders and would in fact be the 
most sensitive one to use in that case. 


B. CP conserving couplings 


The CP conserving dipole couplings interfere with the SM but this contribution 
to the cross-section is suppressed by the r lepton mass as well. It is also possible 
to hnd terms that are linear in the CP conserving anomalous couplings and that 
are not helicity suppressed by looking at single spin correlations. One such term is 
given by 


Oupin ^ af Pa {s{i - u) (pi - P 2 ) ■ (s^- + S^+) + {i - u)^ (Pr+Sp- - Pr-5r+)) 

( 16 ) 

To study Eq. Hn] using only the beam and muon momenta we found the following 
two observables: the hrst one is the muon charge asymmetry 33| dehned by 


Oc = ^\y\ = \yfp+\-\yf,-\ (17) 

The charge asymmetry is C-odd and therefore changes sign at pp collider and van¬ 
ishes at a pp collider as seen in the example in Table IVIHl The second possibility is 
simply 


^Pt qbeam ' {,Pp+ Pp~ ) qbeam ' i^Pp+ T Pp.~ ) 


( 18 ) 


which can also be written as the difference in transverse momentum of the two 
muons. 

To measure any of the correlations discussed above we use the fully integrated 
counting asymmetries normalized to the standard model cross-section. 


A, 


[n+ + n_J 



(19) 
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where iV_|_ = a{Oi > 0), and N_ = a{Oi < 0). The normalization to the standard 
model cross-section, instead of the total cross-section inclnding new physics, is cho¬ 
sen becanse it shows clearly whether a given asymmetry is linear or qnadratic in 
the new conplings. Of conrse, if there is new physics large enongh to be detected 
from deviations in the cross-section from its standard model valne it is possible to 
simply scale the asymmetries we calculate here. It is also possible to look closer at 
details of the angular distributions to attempt to extract new physics, but here we 
limit ourselves to the overall counting asymmetries. 


III. NUMERICAL STUDY 


For our numerical study we generate multiple event samples for the process 
pp —)■ T~^T~ —)■ (where i = fi,e but will be a muon for most of our 

study) at 14 TeV center of mass energy that we summarize in the Appendix. The 


anomalous couplings are implemented in MadGraphB 3^ with the aid of FeynRules 
We use the resulting UFO model hies to generate events for several values of 
drw, drB, drG and d^Q in a range motivated by our previous results from Ref. [13 


The events preserve all spin correlations between production and decay of the r- 
leptons as they are generated for the complete process. In each case we generate 
event samples with one million dimuon or dilepton events after cuts, implying a la 
statistical sensitivity to all asymmetries at the [aja sm x 0.1)% level. 


A. High energy dilepton pairs 

The cuts used in our event generation are: 

• nirT > 120 GeV implemented in the cuts.f hie. The purpose of this cut is to 

exclude the Z resonance region from consideration, as this will be discussed 
separately. We use this idealized cut for simplicity although it may not be pos¬ 
sible to implement experimentally for leptonic tan decays. In a more realistic 
simulation removal of the Z region can be ehectively accomplished with an 
alternative cut on A few tests suggested the cuts give similar results 

but the more realistic one requires much longer event generation time. 

• PTi > 15 GeV for both muons and electrons. This is a standard acceptance 
cut in the LHG experiments. The asymmetries due to new physics increase 
with an increasing pTe cut at the cost of statistical sensitivity. The number 
we use is a good compromise for million event samples. 

• \pe\ < 2.4 is a standard acceptance cut for muons and electrons at the LHG 
experiments and we have checked that this choice does not signihcantly affect 
our asymmetries. 


^ The code is available from the authors upon request. 
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In order to preserve the spin correlations it is important to calculate matrix 
elements for the full process. The numerical implementation of this calculation is 
signihcantly complicated by the very narrow r-lepton width. It is however possible 
to do it with the current version of Madgraph, which has the necessary numerical 
precision, at the cost of long event generation times. A simple trick to alleviate this 
problem is to use a hctitious (and much larger) r-lepton width during the event 
generation, and to then rescale the resulting cross-sections by the narrow-width 
approximation factor VT—fict/^T for each r propagator explicitly 


a 


(^(Xr-fict) X 


^T—fict 

r 

T—exp 


X 


r._ 


SM 




( 20 ) 


The hrst factor is the cross-section calculated by MadGraphB using as an input Tr-fict, 
typically 2.27 x 10“® GeV. The second factor corrects this by 10^^ by rescaling to 
the experimental width. Finally, the last factor takes into account the dependence 
of the r-width on drw- To check that this trick does not distort the kinematic 
distributions of hnal state leptons to the extent of affecting our asymmetries, we 
repeated a few calculations using different hctitious values of the r width spanning 
several orders of magnitude. We show the results of this exercise in Table |Xl 

With the tables presented in the Appendix, we obtain the following approximate 
numerical hts for the most relevant observables: 


- = 1 + 0.019 0.0018 0.0053 Re(d^VK) 

O'SM 

Ai = —0.014 \m.{dT-w) ~ 0.0021 Im(dT-B) 

A 2 = 0.010 \m{drw) + 0.0023 Im(d,-s) 

Ass = 0.0025 Re(dT-iy) \m.{drw) + 0.031 Re(dT-G) Ini(dT-G) + 0.031 Re(d^g) Im(d.,.(5) 
Ac = -0.1125 - 7.1 X 10"^ Re(d^w) 

Ap^ = -0.0955 - 7.0 X 10"^ Re{drw) (21) 


The salient features of these hts are summarized below. 


As discussed in Ref. 1^, terms in the cross-section linear in the real part of 


the anomalous couplings are suppressed by the r mass at LHC energies and 
this is conhrmed both by our ht and by the symmetry of Figured! Our hts are 
not quite the same as the ones we presented in Ref. 1^ due to the diherent 


pre cut used there and the diherent normalization for drw- With the cuts used 
here, the bounds placed on the anomalous couplings by measurements of the 
cross-section (using the same procedure as in Ref. [l3|) are shown in Figured! 
Taking only one parameter to be non-zero at a time we hnd. 


|Im(d^iy)| < 2.7 
—2.9 Xj Re((iT-w) ~ 2.6 


\lm{drB)\ X 8.8 

|Re(d^B)| ^ 8.8 


^ We thank Olivier Mattelaer for this suggestion. 


( 22 ) 












FIG. 1: Regions of drv (left) and the corresponding dr'^, (right) allowed by a maxi- 
mnm 14% deviation from the SM cross-section with the cuts described in the text. 


or equivalently, 


~ 0.0015 , \mT-d'l\ < 0.004 

-0.0016 < < 0.0018 , -0.0068 < m.,a^ < 0.0076 (23) 

These constraints are similar to those quoted in Table [H which is not surprising 
because in both cases they correspond to assuming that the cross-section will 
be measured to 14% accuracy."^. 

• A glance at the Feynman diagrams for qq t^t~ —?• reveals 

that the cross-section should be quadratic in drs and a polynomial of order 
6 in drw because the latter also appears in the r decay vertex as implied by 
the gauge invariant form of the operators, Eq. [31 Our numerical calculation 
indicates that the cross-section has a sensitivity to d^-w at most quadratic, 
in other words the precision of our simulations makes it difficult to allow 
for the higher order terms. This is because our procedure is a form of the 
narrow width approximation (but keeping spin correlations): the dependence 
of a{pp —)■ T^T~) on drw is quadratic, and the r-lepton branching ratios 
remain approximately independent of d^w- 

Interestingly, the r-width itself depends on drw and we could use that to hnd 
an additional constraint. In the approximation in which we treat the hadronic 
r-decay as decay into free quarks, we hnd 

rr(drw) ~ ^t-sm (1 + 0.00126 Ile{dTw) + • • •) (24) 

where • • • stands for much smaller quadratic corrections, and this is the precise 
factor we use in Eq. [201 Of course this approximation does not calculate the 
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This precision corresponds to the largest systematic error in the CMS analysis of high invariant 
mass r-pairs Q 
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hadronic decay modes correctly, but it gives us an estimate^ for the size of the 
corrections introduced by Taken literally, and using from the particle 


data book that the r mean life is (290.3 ± 0.5) x 10 ^^s [3J], it implies 


|Re(dT-iv)| ^ 1.4 


(25) 


• The T-odd correlations Ai ^2 exhibit a linear dependence on the imaginary part 
of the anomalous couplings that is not suppressed by the r mass, as is expected 
for single spin correlations. We hnd that this process is about six times more 
sensitive to \m.{dT-w) than to \m.{dT-B)- 

• The T-odd and CP-odd asymmetry receives contributions quadratic in 
the anomalous couplings that are not suppressed by the r mass. As discussed 
above, they arise from double spin asymmetries produced in the interference 
of the new physics amplitudes with themselves. 

• The independence on the r-mass for two asymmetries originating in the inter¬ 
ference between new physics \m{drw) = 10 and the SM (since Re{dr w) = 0) 
is shown in Table IVIII 

• The T and CP-even asymmetries Ac,pj, exhibit the linear dependence on the 
real part of the anomalous couplings implied by the single spin correlation. 

• At the level of our study, the hnal dimuon channel can be replaced with 
the dilepton channel (including muons and electrons) and this increases the 
statistics by a factor of four without affecting the asymmetries. We explicitly 
compute the asymmetries for the dilepton channel for one value of drw in 
Table IIXI obtaining the same asymmetries as in the dimuon channel given in 
the other Tables. The asymmetries in the dilepton channel are generalized 
from the dimuon case based on the lepton charge (so we also include the 

and hnal states). 


The statistical sensitivity to any of the asymmetries in the dilepton channel with 
100 fb“^ is 0.005. This translates into the following future constraints 


\lm{drw)\ ~ 0.36 , 


m. 


-d5l < 2 X 10 


-4 


|Re(dT-iv)| ~ 7 , 

|Re((iT-vi/)Ini(dTw) I 


\lm{drB) \ ^ 2.2 
ImT-dll < 1 X 10“^ 

I 0.0043, \mra'l\ ^ 0.019 
|Re«G,G)Im«G,G)l ~ 0.16 


(26) 


It is worth noting here that whether we include Eqj2l] or not in Eqj20] only affects 
the constraint we hnd for 'Re{drw), increasing it to 7 from 5. 


® Recall that this approximation used for the SM results in a r-width that is only about 10% 
smaller than the experimental width. 
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B. Dilepton pairs in the Z-resonance region 


The Z-resonance region is selected with the cut 60 < nirT < 120 GeV, with 
the same caveats as before. Keeping the remaining cuts unchanged and generating 
additional samples we obtain the following approximate fits. 

^ = 1 + 0.0069 Re{drwf + 0.0066 \m{drwY + 0-0006 \d^B? + 0.0024 Re{d^w) 

O'SM 

Ai = — 0.013 \m.{drw) — 0.0038 Im^drs) 

A 2 = 0.0016 \m.{drw) + 0.0005 lm{drB) 

= 0.0011 Re{dT-w) lm{drw) + 0.001 Re((iT-G) Im((i.rG) + 0.001 Re((i^g) lm{d^Q) 
Ac = -0.0279 - 5.4 x 10"^ Re(d^w) 

Ap^ = - 0.0251 - 4.7 X 10 "^ Re{drw) ( 27 ) 

We have not included Re{drB) for the asymmetries due to the smaller sensitivity 
already observed in the previous case. Our main observations in this case are: 


Constraints arising from the cross-section are shown in Figured In this case 
we have assumed the cross-section can be measured to 7% accuracy, the current 
systematic uncertainty, following Ref. 35j. Taking only one parameter to be 


15 - 


B 

W 


0.006 - rrird^ 



15 



FIG. 2: Regions of drv (left) and the corresponding (right) allowed by a maxi¬ 

mum 7% deviation from the SM cross-section with the cuts described in the text. 


non-zero at a time we hnd, 

|Im(d.rn^)| < 3.3 , |Im(d^B)| < 10.7 
-3.5 < Re{drw) ^ 3.0 , |Re(d,B)| < 10.6 (28) 

or equivalently, 

\mrd^\ < 0.0018 , \mrd1\ < 0.005 
-0.0018 < < 0.0021 , -0.0078 < < 0.0093 (29) 

which are between 20-25% weaker than those that can be obtained from high 
energy pairs and at best comparable to the existing constraints from LEP. 
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• Assuming again 100 fb“^ implies a much better statistical sensitivity of 0.0009 
when using both electron and muon channels. This better sensitivity is of 
course due to the much larger cross-section and results in the potential con¬ 
straints 


|lm(ci.rvv)| ^ 0.07 , \lm{drB)\ ^ 0.24 


1 < 4.3 X 10-® , 

\mrd1\ < 1.9 X 10 

fe(<^riv) ~ 1-7 , 

mT-af < 0.0011 , ruT-a!) 


|Re(ci^ty)Im(d^v^)| < 0.82 , |Re(ci^g g)Im(ci^g g)| < 0.9 




(30) 


C. Background 

We end this section with a brief discussion of background and how it would 
affect the constraints estimated so far. For the dilepton channel in r-pair production 
both r-leptons in the pair undergo leptonic decay into muons or electrons: pp —)■ 
T^T~ ^ where the missing transverse energy, is due to invisible 

neutrinos. If one lepton is a muon and the other one an electron the dominant 
background arises from tt, W~^W~ or ZZ production. If the two leptons have the 
same flavor there is an additional direct Drell-Yan production of In addition, 

as discussed for example in Ref. 0], contributions from processes where a jet or a 
photon is misidentihed as a lepton are very small. 

The different handles to control this background have been identihed by the 
experimental collaborations. Requiring a minimum missing Et can effectively re¬ 
move the direct Dell-Yan background, we will use > 20 GeV. To suppress ti 
background experiments require at most one jet and no b tagged jets. These re¬ 
quirements are hard to implement at the level of our analysis but they should be 
kept in mind. The requirement that the two leptons be back to back in the trans¬ 
verse plane provides additional suppression against top-pairs and W and Z pairs. 
This is implemented as [i| 


cos A0(r, £+)<-0.95 (31) 

where , £~^) is the difference in azimuthal angle between lepton pairs. And to 

further suppress the contamination from W products, events are selected with an 
additional requirement that the signature being consistent with that of a particle 
decaying into two r leptons. With the following projection variables 

Pg*® = pti+ ■ i + Pti- ■ i, 

pg = pf + • e (32) 

we require p^ — (1.25 x p^*®) > —10, where is a unit vector along the bisector of 
the momenta of the two leptons. 
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Generating MonteCarlo samples for each of the background processes with Mad- 
graph and applying all the preceding cuts to background and signal samples results 
in cross sections: 


a{pp T~^T~ —)■ pT 

a{pp —)■ —)■ Wp'^ p~ $t) 

a{pp —)■ W~^W~ p'^p~$rp) 

a{pp —)■ ZZ —)■ p'^ p~ $rp) 


20.28 fb 
120.1 fb 
18.16 fb 

1.75 fb (33) 


These numbers allow us to quantify the effect of background as follows. First, the 
cuts needed to isolate the signal reduce its cross-section by a factor of 4.7 which 
results in a factor of 2.2 loss in statistical sensitivity. In addition the hnal event 
sample will contain background events that, under our assumptions, are not affected 
by the new physics. If we use the CMS |36| b-tagging efficiency between 70-85% we 


expect less than 11 fb of a{pp —)• tt) background to remain. In this case the T-odd 
asymmetries are reduced by about 2.5. 


Bounds estimated from the cross-sections do not depend so much on the back¬ 
ground cross-section as on its uncertainty and this has already been taken into 
account when we use the experimental estimates for the precision they can achieve 
in their cross-section measurements. The effect of background on the T-even asym¬ 
metries is much harder to estimate, but these do not improve the bounds on Re{drv) 
signihcantly over bounds obtained from cross-sections in any case. 


IV. SUMMARY 

We have examined the possible limits that can be placed on certain anomalous 
couplings of r-leptons at the LHC14 with 100 fb“^. We have considered the four 
dipole-type couplings that appear at dimension six in the effective Lagrangian as 
well as the two r-gluon couplings that appear at dimension eight. We have found the 
statistical sensitivity of single and double spin asymmetries in the dilepton channel 
to these couplings and compared them to the statistical sensitivity from measur¬ 
ing deviations from the SM cross-section at the 14% level. We hnd that T odd 
asymmetries can improve the bounds on the CP violating couplings but that single¬ 
spin asymmetries do not seem to improve the bounds on the anomalous magnetic 
moments. 
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Appendix A: Tables 


All the tables are produced from million event samples for the process pp —?• 
r+r“ at 14 TeV obtained with the following cuts: px^ > 15 GeV, 

IhU — > 120 GeV. With the exceptions noted explicitly below, the r width 

was set to 2.27 x 10“® GeV and the resulting cross-sections were then scaled as 
described in the main text. 


Re((irVK) 

lm{drw) 

cj(fb) 

Ai 

A 2 

^test 

Ac 

A 

■^Pt 

0 

0 

95.45 

0.0003 

0.0000 

0.0019 

-0.1125 

-0.0955 

0 

2 

102.7 

-0.0283 

0.0204 

-0.0002 

-0.1133 

-0.0956 

0 

4 

124.3 

-0.0579 

0.0399 

-0.0013 

-0.1119 

-0.0962 

0 

6 

160.4 

-0.0858 

0.0613 

-0.0018 

-0.1134 

-0.0938 

0 

8 

210.9 

-0.1117 

0.0780 

0.0022 

-0.1088 

-0.0965 

0 

10 

276.0 

-0.1457 

0.1019 

0.0013 

-0.1131 

-0.0991 

2 

0 

103.7 

-0.0003 

-0.0013 

0.0009 

-0.1141 

-0.0964 

2 

2 

111.0 

-0.0294 

0.0194 

0.0013 

-0.1170 

-0.0985 

2 

4 

132.6 

-0.0579 

0.0384 

0.0007 

-0.1175 

-0.0988 

2 

6 

168.6 

-0.0871 

0.0622 

-0.0006 

-0.1150 

-0.0964 

2 

8 

219.4 

-0.1145 

0.0797 

-0.0008 

-0.1171 

-0.1026 

2 

10 

284.2 

-0.1453 

0.1003 

-0.0012 

-0.1206 

-0.1029 

4 

0 

126.4 

-0.0026 

0.0011 

-0.0007 

-0.1155 

-0.1000 

4 

2 

133.6 

-0.0284 

0.0212 

-0.0003 

-0.1178 

-0.1019 

4 

4 

155.3 

-0.0565 

0.0380 

0.0003 

-0.1177 

-0.1009 

4 

6 

191.4 

-0.0837 

0.0640 

-0.0014 

-0.1152 

-0.0986 

4 

8 

242.0 

-0.1140 

0.0820 

-0.0022 

-0.1175 

-0.1012 

4 

10 

306.9 

-0.1445 

0.0971 

0.0022 

-0.1186 

-0.1027 


TABLE II: Single spin T-odd correlations Ai ^2 and T-even correlations Ac,pj, for several 
values of Re((iT-w) and \m.{drw)- Atest should vanish in all cases and gives us an estimate 
of the statistical error. 

In Table [ITl we compute the single-spin asymmetries chosen above for a series 
of values of Re{drw) and Imi^drw) along with Atest which should be zero up to 
statistical error. In Figure |3] we plot the T-odd asymmetries which exhibit the 
expected behavior linear in Im((i.rw)- The hgure also suggests that they have very 
small dependence on Re((iT-iy). In Figure 0] we plot the T-even asymmetries, the 
situation is less clear in this case and we need to study other tables to reach any 
conclusions. 

In Table IIIII we set R.e{drw) = 0 and fit the T-odd correlations to a linear 
equation. This is shown in Figure [5] and the hts are consistent with the previous 
ones from Table UTl We tabulate Atest again to assess the size of statistical fluctuations 
using an asymmetry that should be zero. Finally we also tabulate results for the 
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FIG. 3: A linear fit to lTa.{drw) for Ai (left) and A 2 (right) in Table HIl is supported by the 
data. The separation between the three points corresponding to three values of Iie{drw) 
suggests a very small contribution of the form 'Re{drw)^m.{drw)■ 


- 0.105 


^ - 0.115 


MC ' 


- 0.125 


12 3 4 





FIG. 4: Values of Aq (left) and Ap^ (right) in Table HIl The separation between the 
five points corresponding to the values of lm.{dTw) for each Ke{drw) suggests sizeable 
contributions of the form (Im((i,- w)f- 

T-even asymmetries which should not have a linear dependence on \m{drw)- This 
is conhrmed in Figure [6] within our statistical uncertainty. 

Next we repeat the previous exercise but for lm{drB) instead. The T-odd asym¬ 
metries are also linear in this coupling as expected, but smaller than those induced 
by \v[i{drw) as can be seen in Figure [71 Figure [HI illustrates that the T-even asym¬ 
metries are not affected by this coupling within our numerical sensitivity. 

We turn our attention to the CP conserving couplings in Tables El and EB We 
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l-m{drw) 

(T(fb) 

Til 

A 2 

Atest 


4 

/±Pt 

0 

95.45 

0.0003 

0.0000 

0.0019 

-0.1125 

-0.0955 

2 

102.7 

-0.0283 

0.0204 

- 0.0002 

-0.1133 

-0.0956 

4 

124.3 

-0.0579 

0.0399 

-0.0013 

-0.1119 

-0.0962 

6 

160.4 

-0.0858 

0.0613 

-0.0018 

-0.1134 

-0.0938 

8 

210.9 

-0.1117 

0.0780 

0.0022 

-0.1088 

-0.0965 

10 

276.0 

-0.1457 

0.1019 

0.0013 

-0.1131 

-0.0991 

12 

355.4 

-0.1756 

0.1131 

-0.0041 

-0.1086 

-0.0903 

14 

449.2 

-0.1965 

0.1334 

0.0011 

-0.1135 

-0.0923 

16 

557.0 

-0.2213 

0.1693 

-0.0065 

-0.1087 

-0.0924 

18 

680.0 

-0.2641 

0.1745 

0.0002 

-0.1087 

-0.0886 

20 

817.0 

-0.2726 

0.1793 

-0.0046 

-0.1109 

-0.1007 


TABLE III: Single spin T-odd correlations Ai ^2 and T-even correlations Ac^pj. for several 
values of lm.{drw) with Re(dT-w) = 0 . Atest should vanish in all cases and gives us an 
estimate of the statistical error. 




FIG. 5: MC simulation compared to a linear fit to \m[drw) for Ai (left) and A 2 (right) for 
high rurr events from Table HTTl and also for events in the Z-region 60 < rrirr < 120 GeV . 


also find the expected behavior here: Figure |9] shows that the T-even asymmetries are 
linear in Re((i^vK )5 and Table IVl shows that the T-odd asymmetries are all consistent 
with zero. The dependence of the T-even asymmetries on Iie{drB) is not observed 
within our statistical sensitivity as shown in Figure ITOl 

In Table rvTTl we study the dependence of some of the observables on the r-lepton 
mass. In order to keep kinematic factors in the r decay constant, we also set the 
bottom-quark mass to be always higher than the r mass and we take the charm, 
strange, up and down quarks as well as the muon to be massless. We see that 
the Drell-Yan cross section is approximately independent of the r-lepton mass, as 
expected. The width of the r-lepton exhibits the dependence predicted by the 
SM when the muon mass is neglected. The fact that Ai is approximately constant 
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FIG. 6: MC simulation of ^c,pt high nirT events from Table IIIII and also for events in 
the Z-region 60 < m-T-r < 120 GeV as a function of lm{drw)- 


lm{drB) 

cj(fb) 

4li 

A2 

^test 


A 

jri.prp 

0 

95.45 

0.0003 

0.0000 

0.0019 

-0.1125 

-0.0955 

2 

96.16 

-0.0067 

0.0046 

0.0015 

-0.1129 

-0.0973 

4 

98.25 

-0.0083 

0.0086 

0.0018 

-0.1135 

-0.0959 

6 

101.7 

-0.0122 

0.0143 

0.0016 

-0.1129 

-0.0961 

8 

106.5 

-0.0179 

0.0185 

0.0015 

-0.1141 

-0.0964 

10 

112.7 

-0.0207 

0.0228 

0.0007 

-0.1150 

-0.0982 

12 

120.3 

-0.0233 

0.0288 

0.0012 

-0.1130 

-0.0976 

14 

129.3 

-0.0264 

0.0311 

-0.0008 

-0.1145 

-0.0990 

16 

139.7 

-0.0357 

0.0372 

-0.0013 

-0.1122 

-0.0967 

18 

151.3 

-0.0382 

0.0390 

-0.0013 

-0.1154 

-0.0964 

20 

164.5 

-0.0450 

0.0459 

0.0007 

-0.1148 

-0.0979 


TABLE IV: Single spin T-odd correlations Ai ^2 and T-even correlations Ac^p^ for several 
values of lm{drB) with Re(dT-_B) = 0. A^est should vanish in all cases and gives us an 
estimate of the statistical error. 

in this table supports our interpretation of this result as originating mainly in the 
single spin asymmetry. 

In Table rvTTTl we illustrate the CP properties of the T-odd asymmetries discussed 
above. The asymmetry Ai is GP-odd for the case of the pp collider whereas A 2 is 
GP-even and therefore cannot be induced by the anomalous coupling lm{drw)- For 
the LHC, a pp collider, both T-odd asymmetries are possible as they transform into 
asymmetries in a pp collider under a CP transformation. They do so with opposite 
signs as can be seen in the Table. The charge asymmetry is G-odd and therefore 
changes sign at pp collider and vanishes at a pp collider as seen in the example 
in Table IVIIIl Atest on the other hand is not zero for pp colliders where the beam 
direction can be dehned unambiguously. 
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FIG. 7: MC simulation compared to a linear fit to lm{drB) for (left) and A 2 (right) for 
high TUrr events from Table HVl and also for events in the Z-region 60 < m-j-r < 120 GeV . 
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FIG. 8: MC simulation of Ac^pj, for high rrirr events from Table IIVI and also for events in 
the Z-region 60 < uirr < 120 GeV as a function of lm{drB)- 

Table HXl illustrates that replacing the dimuon channel with the dilepton channel 
simply increases the statistics by a factor of four and does not change the four 
asymmetries we have been discussing. 

In Table |X] we demonstrate that the trick of using a hctitious r-lepton width in 
the simulations does not affect the cross-section or the asymmetries A 2 and Ac (it 
also does not affect the other asymmetries). 

In Table IXll we set Re{drw) = R^{drw) to look for the double spin asymmetry 
through The result of the ht is as shown in Figure [TTl In Table Kill we set 
Re(dT-G) = Im((iT-G') to look for the double spin asymmetry through The result 
of the ht is as shown in Figure [TTJ 

In Table IXim we set Re((i^g) = Im((i^( 5 ) to look for the double spin asymmetry 
through Ass. The result of the ht is as shown in Figure [TTl 
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Ke{drw) 

cj(fb) 

4li 

A2 

^test 


A 

i±Pt 

-20 

806.2 

-0.0059 

-0.0033 

-0.0208 

-0.1022 

-0.0899 

-16 

548.8 

-0.0045 

0.0007 

0.0034 

-0.0998 

-0.0876 

-12 

349.3 

0.0000 

-0.0017 

-0.0034 

-0.1008 

-0.0807 

-8 

206.9 

0.0015 

0.0013 

0.0004 

-0.1078 

-0.0915 

-4 

122.3 

0.0017 

-0.0001 

0.0002 

-0.1105 

-0.0935 

0 

95.45 

0.0003 

-0.0000 

0.0019 

-0.1125 

-0.0955 

4 

126.4 

-0.0026 

0.0011 

-0.0007 

-0.1155 

-0.0999 

8 

215.1 

-0.0020 

-0.0005 

0.0001 

-0.1235 

-0.1060 

12 

362.1 

0.0045 

0.0000 

0.0014 

-0.1236 

-0.1023 

16 

566.9 

0.0039 

0.0036 

0.0017 

-0.1273 

-0.1105 

20 

830.3 

0.0145 

0.0079 

0.0019 

-0.1221 

-0.1127 


TABLE V: Single spin T-odd correlations Ai ^2 and T-even correlations Ac,pj. for several 
values of Re{drw) with Im(d^vu) = 0- Atest should vanish in all cases and gives us an 
estimate of the statistical error. 






FIG. 9: MC simulation compared to a linear fit to Re(dT-w) for Ac (left) and Ap^, (right) 
for high nirT events from Table IVl and also for events in the Z-region 60 < nirT < 120 GeV. 
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(j(fb) 

Ai 

A 2 

Afest 

Ac 

A 

^Pt 

-20 

165.9 

-0.0002 

-0.0007 

0.0016 

-0.1155 

-0.0950 

-16 

140.7 

-0.0007 

-0.0028 

-0.0006 

-0.1139 

-0.0981 

-12 

121.1 

0.0001 

0.0003 

-0.0017 

-0.1142 

-0.0962 

-8 

107.1 

0.0008 

-0.0010 

0.0005 

-0.1112 

-0.0950 

-4 

98.54 

0.0005 

0.0000 

0.0005 

-0.1128 

-0.0956 

0 

95.45 

0.0003 

-0.0000 

0.0019 

-0.1125 

-0.0955 

4 

97.95 

-0.0000 

-0.0016 

-0.0008 

-0.1126 

-0.0958 

8 

105.9 

-0.0009 

-0.0000 

0.0014 

-0.1135 

-0.0962 

12 

119.5 

0.0003 

0.0002 

-0.0002 

-0.1137 

-0.0983 

16 

138.6 

0.0017 

-0.0004 

0.0002 

-0.1147 

-0.0981 

20 

163.1 

-0.0035 

-0.0026 

-0.0026 

-0.1176 

-0.1002 


TABLE VI: Single spin T-odd correlations Ai ^2 and T-even correlations Ac,pj. for several 
values of Re(dT-B) with Im(dT-B) = 0. Atest should vanish in all cases and gives us an 
estimate of the statistical error. 
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-0.113 
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FIG. 10: MC simulation compared to a quadratic fit to Re(dT-B) for Aq (left) and Ap^ 
(right) in Table IVll 
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mr(GeV) 


Rr(r+ 

a{T~^T ) 

v's) 



wSaSmM 

2.5 

- 

0.123 

21.53 pb 

138.9 fb 

-0.0721 

-0.0589 

3.0 

1.002 

0.124 

21.53 pb 

139.0 fb 

-0.0696 

-0.0561 

3.5 

1.001 

0.125 

21.53 pb 

139.7 fb 

-0.0709 

-0.0576 

4.0 

1.002 

0.125 

21.52 pb 

139.6 fb 

-0.0729 

-0.0549 

4.5 

1.002 

0.125 

21.50 pb 

139.7 fb 

-0.0725 

-0.0585 

5.0 

1.003 

0.125 

21.48 pb 

139.6 fb 

-0.0697 

-0.0573 

5.5 

1.004 

0.125 

21.46 pb 

139.5 fb 

-0.0710 

-0.0574 

6.0 

1.005 

0.125 

21.45 pb 

139.4 fb 

-0.0724 

-0.0600 

6.5 

1.006 

0.125 

21.43 pb 

139.2 fb 

-0.0736 

-0.0564 

7.0 

1.007 

0.125 

21.40 pb 

139.2 fb 

-0.0727 

-0.0583 


TABLE VII; Scaling with the r mass of cross-section and two asymmetries Ai and Ac 
for Im((iT-vy)=10. To remove any kinematic dependence from the decay vertices we use 
niu = rud = rric = rus = 0, nib = 10 in all cases. 


Collider 

cj(fb) 

Til 

A 2 

^test 


pp 

276.0 

-0.1457 

0.1019 

0.0013 

-0.1131 

pp 

275.8 

-0.1418 

-0.0999 

-0.0021 

0.1177 

pp 

313.6 

-0.1531 

0.0021 

0.1687 

-0.0005 


TABLE VIIL Comparison of T-odd and T-even asymmetries with Re(dT-w)=0, 
Im((iT-vy) = I0 for different colliders to exhibit their transformation properties under CP. 


lm{drw) 

cj(fb) 


A 2 


A 

■^Pt 

10 

1111.0 

-0.1363 

0.0974 

-0.1100 

-0.0956 


TABLE IX: Selected asymmetries in the dilepton channel {pp —)■ r^r ^ n's) with 
Re((iT-iv)=0. In this case itsm = 385.2 fb. 


r^(GeV) 

lTa.{drw) 

cj(fb) 

A 2 

Tic 

2.27 X 10“^ 

10 

276.0 

0.1019 

-0.1131 

2.27 X 10"® 

10 

276.0 

0.0952 

-0.1156 

2.27 X 10“^ 

10 

275.8 

0.0960 

-0.1179 

2.27 X 10“® 

10 

275.9 

0.0953 

-0.1169 

2.27 X 10-*^ 

10 

276.0 

0.0975 

-0.1155 


TABLE X; Effect of changing the r-lepton width in the MC simulation. After the rescaling 
described in the text the cross-section as well as the asymmetries are seen to be indepen¬ 
dent of Tt within our numerical accuracy. In all cases we took Re((iT-vu)=0. 
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Re((iT-vu) 

l'ni{drw) 

cj(fb) 


3 

3 

129.5 

0.0209 

4 

4 

155.3 

0.0399 

5 

5 

188.4 

0.0623 

6 

6 

228.6 

0.0919 

7 

7 

276.2 

0.1177 

8 

8 

331.0 

0.1544 

9 

9 

392.9 

0.1959 

10 

10 

462.2 

0.2483 

11 

11 

538.5 

0.3095 

12 

12 

622.2 

0.3660 


TABLE XI; Double spin correlation induced by interference between Jie{drw) and 
lm.{drw)- 


Reidrc) 

Im((i^G) 

cj(fb) 



A 

■^Pt 

0.4 

0.4 

105.1 

0.0039 

-0.1123 

-0.0955 

0.8 

0.8 

134.1 

0.0204 

-0.1152 

-0.0975 

1.2 

1.2 

182.2 

-0.0455 

-0.1126 

-0.0966 

1.6 

1.6 

249.8 

0.0791 

-0.1118 

-0.0948 

2.0 

2.0 

336.7 

0.1232 

-0.1084 

-0.0923 


TABLE XII: Double spin correlation induced by interference between Iie{drG) and 
Im((iT-G'). No discernible effect from these couplings is found in other asymmetries. 


Re(d^^) 

Ini«G) 

cj(fb) 



A 

■^Pt 

0.4 

0.4 

105.1 

0.0049 

-0.1132 

-0.0945 

0.8 

0.8 

134.1 

0.0193 

-0.1131 

-0.0964 

1.2 

1.2 

182.2 

0.0455 

-0.1124 

-0.0974 

1.6 

1.6 

249.8 

0.0778 

-0.1155 

-0.0971 

2.0 

2.0 

336.7 

0.1234 

-0.1117 

-0.0951 


TABLE XIIL Double spin correlation A^^ induced by interference between Re((i^g) and 
Im((i^^). No discernible effect from these couplings is found in other asymmetries. 
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FIG. 11: Comparison of MC and a fit linear in Yle{drw)^"^{dTw) (up), Re((iT-G)Iui((iT-G), 
(left) and Re((i.^^)Im((i^g) (right), for shown separately for events with high nirT and 
for events in the Z-region 60 < ttIt-t < 120 GeV . 


Re(drG) 

1—1 

cj(fb) 


^G 

A 

i±Pt 

2.0 

2.0 

336.9 

-0.0058 

-0.1116 

-0.0956 

Re((i^(5) 

1—1 

cj(fb) 


^G 

A 

^Pt 

2.0 

2.0 

336.9 

-0.0026 

-0.1151 

-0.0987 


TABLE XIV: The real and imaginary parts of the different couplings drc and d^Q do not 
show interference effects. 
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